(12) INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (POT) 



(19) World Intellectual Property Organization 
International Bureau 

(43) International Publication Date 
29 August 2002 (29.08.2002) 




PCT 



lilllllllllllllliillllili 

(10) International Publication Number 

wo 02/067182 A2 



(51) International Patent Classification 



7. 



G06F 19/00 



(21) International AppUcation Number: PCT/US02/05553 

(22) International Filing Date: 20 February 2002(20.02.2002) 



(25) Filing Language: 



English 



(81) Designated States (national): AE, AG, AL. AM, AT, AU, 
AZ, BA, BB, BG, BR, BY, BZ, OA, CH, CN, CO, CR, CU, 
CZ, DE, DK, DM, DZ, EC, EE, ES, H, GB, GD, GE, GH, 
GM, MR, HU, ID, IL, IN, IS, JP, KE, KG, KP, KR, KZ, LC, 
LK, LR, LS, LT, LU, LV, MA. MD, MG, MK, MN, MW, 
MX, MZ, NO, NZ, PH, PL, PT, RO, RU, SD. SE, SG, SI, 
SK, SL, TJ, TM, TR, TT, TZ, UA, UG, US, UZ, VN, YU, 
ZA, ZW. 



(26) Publication language: 



English 



(30) Priority Data: 

09/789,595 



20 February 2001 (20.02.2001) US 



(71) Applicant (for all designated States except US): CY- 
TOKINETICS, INC. (US/US]; 280 E. Grand Avenue, 
Suite 2, South San Francisco, CA 94080 (US). 

(72) Inventors; and 

(75) Inventors/Applicants (for US only): ADAMS, Cynthia, 
L. fUS/US]; 2409 Cedar Street, Berkeley, CA 94708 (US). 
OESTREICHER, Donald, R. [USAJS]; 904 Old Town 
Court, Cupertino, CA 95014-4024 (US). VAISBERG, 
Eugeni, A. [USAJS]; 647 Pegasus Lane, Foster City, CA 
94404 (US). 

(74) Agent: WEAVER, Jeffrey, K.; Beyer Weaver & Thomas, 
LLP, P.O. Box 778, Berkeley, CA 94704-0778 (US). 



(84) Designated States (regional): ARIPO patent (GH, GM, 
KE, LS, MW, MZ, SD, SL, SZ, TZ, UG. ZM, ZW), 
Eurasian patent (AM, AZ, BY, KG, KZ, MD, RU, TJ, TM), 
European patent (AT, BE, CH, CY, DE, DK, ES, FX, FR, 
GB, GR, IE, IT, LU, MC, NL, PT, SE, TR), OAPI patent 
(BF, BJ, CF, CG, CI, CM, GA, GN, GQ, GW, ML, MR, 
NE, SN, TO, TG). 

Published: 

— - without international search report and to be republished 
upon receipt of that report 

For two-letter codes and other abbreviations, refer to the **Guid' 
ance Notes on Codes and Abbreviations*' appearing at the begin- 
ning of each regular issue of the PCT Gazette. 



(54) Title: CHARACTERIZING BIOLOGICAL STIMULI BY RESPONSE CURVES 



103 



Ids 



IDEhmFY COLLECTION OP 
CHEMCAL CXIMKIUKDS FOR 
USE WAKALV51S 



101 



SEl£CTCUMtEOTCOMK>UND - 



107- 



SET fflEW COMBIMATKJN OF 

DOSE AND CELL UNE 



109 



*\ IMA<3 CELL* OF CURREMT CELL 
^ LINE AT CWWENT DOSE 



ni 



00 



MEASURE AND STORE VALUES 

Of MULTITLE PHENQFTirWC 
FEATURES {AT MUtTVlE TIMES 
IFNECESSADY) 



111 



A^'''l3oRECOMBIKATHlNS OF DoS^vJ 
^•^.^^ AND CELL LtME? 



MS 



NO 



COMBINE FEATURE VALUES ACROSS 

MULTIPLE CELL LDffiS TO OerTAlN 
SEPARATE mENOrnmC VECTORS FOR 
EAOIDOSE 



U1 



o 



119 



^ 



YES 



(57) Abstract: A method for generating stimulus response curves (e.g., dose response 
curves) shows how the phenotype of one or more cells change in response to varying lev- 
els of the stimulus. Each "point" on the curve represents quantitative phenotype forcell(s) 
at a particular level of stimulus (e.g., dose of a therapeutic). The quantitative phenotypes 
are multivariate phenotypic representations of the cell(s). They include various features 
of the cell(s) obtained by image analysis. Such features often include basic parameters 
obtained from images (e.g., cell shape, nucleus area, Golgi texture) and/or biological 
characterizations derived firom the basic parameters (e.g., cell cycle state, mitotic index, 
etc.). The stimulus response curves may be compared to allow classification of stimuli 
and identify subtle differences in related stimuli. To facilitate the comparison, it may be 
desirable to present the response curves in a principal component space. 



REIJUCE WMENSKNALrrV or FEATURE SF ACE 
AND DenCT RESPONSE PAIHS FOR EACH 
DOMKXJND 



121 



COMPARE SEPARATE PATHS OF 
TMECOMPOUKDS 



wo 02/067182 



PCT/US02/05553 



PATENT APPLICATION 

5 CHARACTERIZING BIOLOGICAL STIMULI BY 

RESPONSE CURVES 

BACKGROUND OF THE INVENTION 

The present invention relates to techniques for determining the response of 
10 biological cells to varying levels of a particular stimulus. More specifically, the 
invention relates to response curves derived fi-om multivariate phenotypic data 
extracted from images of biological cells. 

Piffified substances having a desirable combination bio-active properties are 
rare and often difficult to identify. Recent advances in traditional organic chemistry 

15 and the development of rapid combinatorial chemistry techniques have increased the 
number of compounds that researchers can test for a specific biological activity (e.g., 
binding to a target). Unfortunately, the vast majority of "hits" generated by such 
techniques do not possess the right combination of properties to qualify as therapeutic 
compounds. When these substances are subjected to low throughput cellular and 

20 animal tests to establish their therapeutic usefuhiess, they are typically found to fail in 
some regard. Unfortunately, such tests are time consuming and costly, thus hmiting 
the number of substances that can be tested. In a like regard, the few hits that do 
possess the right combination of properties avoid recognition until after the 
throughput tests are conducted. With better early evaluation techniques, such 

25 promising candidates could be identified earUer in the development process and put 
on a fast track to the marketplace. 

Various early evaluation techniques are under investigation and some have 
shown promise. In particular cellular phenotyping technologies employing 
sophisticated image analysis have proven very useful in characterizing thorapeutic 
30 chemicals. Such technologies are generally described in WO/00/70528 published on 
November 23, 2000. These techniques attempt to classify compounds based on 
phenotypic changes that they induce. From these changes, detailed mechanisms of 
action can be deduced. 

Typically, researchers attempting to classify a new compound based on 
35 mechanism of action wish to know how that compound compares to other known 
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therapeutics. Compounds that exhibit similar biological functioning m some regards 
may exhibit similarity in other regards as well. One difficulty in assessing similarity 
is that compoimds often have greatly varying potencies. In other words, while two 
different compounds may operate by the same or similar mechanism of action, one 
5 compound may operate at a much lower concentration than the other compound. It is 
difficult to make meaningful comparison of two such compounds until the dose scales 
of these compounds have been adjusted. To this end, researchers often use dose 
response curves to compare compounds. These curves show the biological 
effectiveness of particular drugs over multiple concentrations. The effect of the dmg 
10 at each different concentration provides the ^'points" for the dose response curves. 

Typically, such dose response curves are limited to a single particular 
biological parameter (e.g., cell count or expression of a protein). The numeric value 
of such parameter is provided as a function of concentration for each compound of 
interest. The resulting curves can be compared to identify similar trajectories. Two 
15 compounds having similar trajectories might be expected to operate by the same 
mechanism of action, depending upon which biological parameter is being considered. 
Unfortunately, there are significant limits to the value of such comparisons. Most 
importantly, many different parameters may contribute to a mechanism's signature. 
So a simple dose response curve may fail to shed Ught on a mechanism. 

20 While image analysis techniques for characterizing phenotypes can provide 

many different characteristics of a compound, their full potential has not yet been 
realized. Particularly, it would be usefiil if such techniques could be applied to obtain 
meaningful dose response information for compounds or other stimuli under 
investigation. 

25 

STTMM ARY OF THE INVENTION 

The present invention provides a method, program code, and apparatus for 
generating stimulus response curves (e.g., dose response curves) showing how the 
phenotype of one or more cells change in response to varying levels of the stimulus, 

30 Each "point" on the curve represents quantitative phenotype for cell(s) at a particular 
level of stimulus (e.g., dose of a therapeutic). The quantitative phenotypes are 
multivariate phenotypic representations of the cell(s). They include various features 
of the cell(s) obtained by image analysis. Such featm-es often include basic 
parameters obtained from images (e.g., cell shape, nucleus area, Golgi texture) and/or 

35 biological characterizations derived from the basic parameters (e.g., cell cycle state, 
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mitotic index, etc.). The stimulus response curves may be compared to allow 
classification of stimuli and identify subtle differences in related stimuli. To facilitate 
the comparison, it may be desirable to present the response curves in a principal 
component space. 

5 One specific aspect of the invention provides a method for determining the 

response of cells to multiple levels of a stimulus. The method may be characterized 
by the following sequence: (a) obtaining feature values which characterize the 
phenotype of cells exposed to a particular level of the stimulus to produce a 
"quantitative phenotype," (b) repeating (a) for each of the multiple levels of stunulus 

10 to thereby produce a separate quantitative phenotype of the cells at each level of 
stimulus; and (c) identifying a path through the separate quantitative phenotypes of 
cells exposed to the stimulus. The stimulus can take many different forms. Examples 
include exposure to chemical compounds, exposure to biological agents, exposure to 
electromagnetic radiation, exposure to particle radiation, exposure to an electrical or 

15 magnetic field or force, exposure to a mechanical field or force, and combinations of 
these. In some cases, the multiple levels of stimulus are multiple durations of after an 
initial exposure to the stimulus. In this embodiment, the cells are analyzed at various 
times after exposure. 

In particularly preferred embodiments, at least some feature values comprising 
20 the quantitative phenotypes are obtained firom an image of the cells. These feature 
values may characterize cell morphology, statistical features of cells (sometunes 
derived fi-om intensity histograms), biological classification of the cells, and the like. 
In one example, a biological classification specifies a cell cycle state. 

Sometimes a graphical representation of the path provides the most usefiil 
25 information. In a particularly prefenred embodhnent, the graphical representation is 
provided along one or more principle components obtained via a principle component 
analysis. 

Another aspect of the mvention pertains to apparatus for analyzing images of 
cells exposed to multiples of a stimulus and generating a response path based on those 

30 images. The apparatus includes at least (a) an interface configured to receive the 
images of the cells that have been exposed to said multiple levels of a stimulus; (b) a 
memory for storing, at least temporarily, some or all of the images; and (c) one or 
more processors in communication with the memory and designed or configured to 
generate the response path by a technique of the type described herein. Typically, the 

35 apparatus will also include a display that is capable of graphically depicting the 
response path. 

3 
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As indicated, the invention provides particular value when used to determine 
whether a first compound and a second compound act on cells by a related mechanism 
of action. Thus, another aspect of the invention may be characterized by the following 
sequence: (a) for each of multiple concentrations of the first compound, obtaining a 
5 plurahty of feature values characterizing the phenotype of cells exposed to the 
particular concentration of the first compound, to thereby produce a plurality of first 
concentration-specific phenotypes; (b) identifying a first path through the first 
concentration-specific phenotypes of cells exposed to the first compound; (c) for each 
of multiple concentrations of the second compound, obtaining a plurality of feature 

10 values characterizing the phenotype of cells exposed to the particular concentration of 
the second compound, to thereby produce a plurahty of second concentration-specific 
phenotypes; (d) identifying a second path through second concentration-specific 
phenotypes of cells exposed to the second compound; and (e) comparing the first and 
second paths, wherein a degree of similarity between the paths corresponds to a 

15 degree of similarity in the mechanism of action of the first and second compounds. In 
some particularly valuable apphcations, at least one of the first and second compounds 
is a known therapeutic or potential therapeutic. 

The concentrations of the compounds should vary over an active range. The 
multiple concentrations of the first compoxmd typically vary firom lowest to highest by 
20 a factor of at least about two. Preferably, the multiple concentrations of the first 
compoimd include at least five separate concentrations of the first compound, and 
more preferably at least eight separate concentrations of the first compound. 

As mentioned above in the context of the first aspect of the invention, the 
feature values may be provided firom a number of different sources. Particularly 

25 valuable phenotypic features are provided by image analysis and associated processes, 
hi a particularly preferred embodiment of this aspect of the invention, the feature 
values include nimieric values characterizing one or more of the following cellular 
components: DNA, Golgi, cytoskeletal components such as tubulin and actin, and the 
plasma membrane. In a specific embodiment, the plurahty of feature values include 

30 numeric values characterizing one or more of the following cellular components: 
DNA, Golgi, and tubulin. 

In one approach, the comparison simply involves graphically depicting the first 
and second paths together. Preferably, the graphical depiction presents the first and 
second paths in a space defined by principal components. Thus, in some 
35 embodiments, the method also involves using the concentration specific phenotypes in 
a technique that provides a reduced-dimensionahty space in which to depict the paths 
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(e.g., principal component analysis, linear and non-linear discriminant analysis, 
multidimensional scaling, and projection pursuit techniques). 

Another aspect of the invention pertains to computer program products 
including machine-readable media on which are stored program instructions for 
5 implementing at least some portion of the methods described above. Any of the 
methods of this invention may be represented, in whole or in part, as program 
instructions that can be provided on such computer readable media. In addition, the 
invention pertains to various combiaations of data and data structures generated 
and/or used as described herein. 

10 These and other features and advantages of the present invention wiU be 

described in more detail below with reference to the associated figures. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a process flow chart depicting the preparation and use of a stimulus 
1 5 response curve based upon phenotypic data. 

Figure 2 is a simplified block diagram of a computer system that may be used 
to ioaplement various aspects of this invention such as liie various image analysis 
algorithms of this invention. 

Figure 3A is a plot of several dose response curves for compounds having 
20 known effects on targets; the ciirves are presented in a space defined by three principal 
components defined for multivariate phenotypic information. 

Figure 3B is a plot of the dose response curves of only the faraesyltransfarase 
and geranylgeranyltransferase inhibitors from Figure 3 A. 

Figure 3C is a plot of the dose response curves of only the actin and tubuUn 
25 inhibitors from Figure 3 A. 

Figure 3D is a plot of the dose response cinrves of only the compoimds that 
effect signaling pathways directly from Figure 3 A. 

Figure 4 is a PCA plot that highlights the deviation of the actin inhibitor 
Cytochalasin A from the other actin inhibitors. 
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Figures 5A and 5B are plots of simple dose response curves of a tubulin 
feature as it varies with concentration of Cytochalasin A and Cytochalasin J, 
respectively. 

Figure 6 is a dendrogram showing the compounds of Figure 3 A at the IC50 for 
5 A549 cells. 

Figures 7 A and 7B are plots of simple dose response curves, across cell lines, 
using cell comt as indicator of the potency of Mastoparan and its synthetic analog 
MAS7, respectively. 

Figure 8 is a PCA plot showing the dose response paths for 83 different 
10 oncology compounds clustered in groups representative of the mechanism of action. 

Figure 9 is a PCA plot showing compounds that have biochemical activity 
against an oncology target described in the examples. 

Figure 10 is a PCA plot showing a subset of the compounds in Figure 9, which 
subset was identified in a primary screen. 

15 Figure 1 1 is a plot shovmig a region of PCA space that represents the optimal 

profile for inhibition of an oncology target. 

Figure 12 is a zoomed-in figure of the compoimds firom Figure 11 in PCA 

space. 

Figure 13 is a process flow chart showing how a single multivariate dose 
20 response experiment was used to identify and narrow down 57 hits to eight highly 
potent and specific compounds. 



DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

The present invention allows for comparison and visualizing of response 
25 curves in multidimensional space. The response curves may span various levels of a 
stimulus, with each point in the curve representing a different level of the stimulus. 
For example, each point might represent a different concentration or dose of chemical 
compound. Alternatively, each pomt in the curve may represent a diJBFerent time after 
initial exposure to a chemical compound. Importantly, each point in the response 
30 curve contains multivariate information about a cell's or population of cells' response 
to a particular level of the stimulus. Preferably, this multivariate inforaiation contains 
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some phenotypic infonnation about the cell. Such phenotypic informatioii may 
provide morphological details, statistical details, and/or higher level biological 
characterizatioiis of the cell or cell population. In an especially preferred 
embodiment, such features are extracted directly or mdirectly from images of the 
5 cells. Of course, the multivariate information in the data points may include non- 
phenotypic information as well. Such information can derive from any of a ntmiber of 
different tests and/or other sources such as pubUc literature and databases. 

One important advantage of the present invention is that it allows related 
stimuli to be compared in a manner that accounts for complicated interactions 

10 between multiple phenotypic variables. Such comparisons help identify trends and 
allow characterization of particular stimuli. The comparisons may be accomphshed 
by a computing device and/or human observers. To the extent that human observers 
are involved in the comparison, it will be beneficial to depict the multivariate response 
curves in a space that emphasizes variations in the data. For example, the invention 

1 5 may involve depicting the response paths in a space defined by principle components. 
In this manner, complicated multivariate data is depicted so that it can be easily 
comprehended. To the extent that a quantitative comparison is required, a computing 
device may compare two or more response curves by any of a number or techniques. 
Such techniques include distance techniques, clustering techniques and the like, 

20 

Process Overview and Relevant Definitions 

Some of terms used herein are not conunonly used in the art. Other terms may 
have multiple meanings in the art. Therefore, the following definitions are provided 
as an aid to understanding the description that follows. The invention as set forth in 
25 the claims should not necessarily be limited by these definitions. 

The term "component" or "component of a cell" refers to a part of a cell 
having some interesting property that can be employed to derive biologically relevant 
information using image analysis. General examples of cell components include 
biomolecules and cellular organelles. Specific examples of biomolecules that could 

30 serve as cell components for use with this invention include proteins, lipids, 
polysaccharides, proteins, etc. Sometimes, the relevant component will refer to a 
group of structurally or fimctionally related biomolecules. Alternatively, the 
component may represent a portion of a biomolecule such as a polysaccharide group 
on a protein, or a particular sequence of a nucleic acid or protein. Collections of 

35 molecules such as micells can also serve as cellular components for use with this 

7 
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invention. And subcellular structures such as vesicles and organelles may also serve 
the purpose. 

The term "marker" or "labeling agent* ' refers to materials that specifically bind 
to and label cell components. These markers or labeling agents shoxild be detectable 
5 in an image of the relevant cells. Typically, a labeUng agent emits a signal whose 
intensity is related to the concentration of the cell component to which the agent 
binds. Preferably, the signal intensity is directly proportional to the concentration of 
the underlying cell component. The location of the signal source (i.e., the position of 
the marker) should be detectable in an image of the relevant cells. 

10 Preferably, the chosen marker binds iadiscriminately with its corresponding 

cellular component, regardless of location within the cell. Allhougb in other 
embodiments, the chosen marker may bind to specific subsets of the component of 
interest (e.g., it binds only to sequences of DNA or regions of a chromosome). The 
marker should provide a strong contrast to other features in a given image. To this 

15 end, the marker should be luminescent, radioactive, fluorescent, etc. Various stains 
and compounds may serve this purpose. Examples of such conqwunds include 
fluorescently labeled antibodies to the cellular component of interest, fluorescent 
intercalators, and fluorescent lectms. The antibodies may be fluorescently labeled 
eithCT directly or iudirectly. 

20 The term "stimulus" refers to something that may influence the biological 

condition of a cell. Often the term will be synonymous with "agent" or 
"manipulation." Stimuli may be materials, radiation (including all maimer of 
electromagnetic and particle radiation), forces (including mechanical (e.g., 
gravitational), electrical, magnetic, and. nuclear), fields, thermal energy, and the like, 

25 General examples of materials tiiiat may be used as stimuli include organic and 
inorganic chemical compovmds, biological materials such as nucleic acids, 
carbohydrates, proteins and peptides, hpids, various infectious agents, mixtures of the 
foregoing, and the like. Other general examples of stimuli include non-ambient 
temperature, non-ambient pressure, acoustic energy, electromagnetic radiation of all 

30 frequencies, the lack of a particular material (e.g„ the lack of oxygen as in ischemia), 
temporal factors, etc. 

Specific examples of biological stimuli include exposure to hormones, growth 
factors, antibodies, or extracellular matrix components. Or exposure to biologies such 
as infective materials such as viruses that may be naturally occurring viruses or 
35 vimses engineered to express exogenous genes at various levels. Biological stimuU 
could also include delivery of antisense polynucleotides by means such as gene 

8 
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transfection. Stimuli also could include exposiu-e of cells to conditions that promote 
cell fusion. Specific physical stimuU could include exposing cells to shear stress 
under different rates of fluid flow, exposmre of cells to different temperatures, 
exposure of cells to vacumn or positive pressure, or exposure of cells to sonication, 
5 Another stimxilus includes applying centrifugal force. Still other specific stimuli 
include changes in gravitational force, including sub-gravitation, appHcation of a 
constant or pulsed electrical current Still other stimuli include photobleaching, which 
in some embodiments may include prior addition of a substance that would 
specifically mark areas to be photobleached by subsequent light exposure. In 
10 addition, these types of stimuli may be varied as to time of exposure, or cells could be 
subjected to multiple stimuli in various combinations and orders of addition. Of 
course, the type of manipulation used depends upon the application. 

The term "phenotype" generally refers to the total appearance of an organism 
or cell from an organism. In the context of this invention, cellvdar phenotypes and 

15 their representations in processing systems (e.g., computers) are particularly 
interesting. A given cell's phenotype is a function of its genetic constitution and 
environment. Often a particular phenotype can be correlated or associated with a 
particular biological condition or mechanism of action resulting from exposure to a 
stimulus. Generally, cells undergoing a change in biological conditions will undergo a 

20 corresponding change in phenotype. Thus, cellular phenotypic data and 
characterizations may be exploited to deduce mechanisms of action and other aspects 
of cellular responses to various stimuli. 

A selected collection of data and characterizations that represent a phenotype 
of a given cell or group of cells is sometimes referred to as a "quantitative cellular 

25 phenotype." This combination is also sometimes referred to as a phenotypic 
fingerprint or just "fingerprint." The multiple cellular attributes or features of the 
quantitative phenotype can be collectively stored and/or indexed, nxmierically or 
otherwise. The attributes are typically quantified in the context of specific cellular 
components or markers. Measured attributes useful for characterizing an associated 

30 phenotype include morphological descriptors (e.g., size, shape, and/or location of the 
organelle) and composition (e.g., concentration distribution of particular biomolecules 
within the organelle). Other attributes include changes in a migration pattern, a 
growth rate, cord formation, an extracellular matrix deposition, and even cell count. 

The quantitative phenotypes may themselves serve as individual points on 
35 response curves of this invention. A phenotypic response to stimulus may be 
characterized by exposing various cell lines to a stimulus of interest at various levels 
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(e.g., doses of radiation or concentrations of a compound). In each level within this 
range, the phenotypic descriptors of interest are measured to generate quantitative 
phenotypes associated with levels of stimulus. 

The tenn "path" or "response curve" refers to the characterization of a 
5 stimulus at varioxis levels. For example, the path may characterize the effect of a 
chemical appUed at various concentrations or the effect of electromagnetic radiation 
provided to cells at various levels of intensity or the eflfect of depriving a cell of 
various levels of a nutrient. Mathematically, the patii is made up of multiple points, 
each at a different level of the stimulus. In accordance with this invention, each of 
10 these points is preferably a collection of parameters or characterizations describing 
some aspect of a cell or collection of cells. Typically, at least some of these 
parameters and/or characterizations are derived from images of the cells. In this 
regard, they represent quantitative phenotypes of the cells. In the sense that each point 
in the path may contain more than one piece of information about a cell, the points 
15 may be viewed as arrays, vectors, matrices, etc. To the extent that the path coimects 
^ points 'containing phenotypic information (separate quantitative phenotypes), the path 
itself may be viewed as a "concentration-independent phenotype." 

As used herein, the term "feature" refers to a phenotypic property of a cell or 
population of cells. Typically, the points in a response curve of this invention are 

20 each comprised of multiple features. The terms "descriptor"' and "attribute" may be 
used synonymously with "feature." Features derived from cell images include both 
the basic "parameters" extracted from a cell image and the "biological 
characterizations" (including biological classifications such as cell cycle states). The 
latter example of a feature is typically obtained from an algorithm that acts on the 

25 basic parameters. The basic parameters are typically morphological, concentration, 
and/or statistical values obtained by analyzing a cell image showing the positions and 
concentrations of one or more markers bomd v^thin the cells. 

Figure 1 depicts a sample process flow for generating and using response paths 
in accordance with an embodiment of this invention. As depicted in Figure 1, a 
30 process 101 begins by identifying a collection of chemical compounds for use in the 
analysis. See block 103. This operation may be performed by a computing apparatus 
or possibly by one or more human beings. The compounds selected at 103 will 
ultimately be used to generate data that defines a "phenotypic space" for comparing 
multiple response paths. 

35 After the relevant collection of chemical compounds has been identified at 

103, the process next selects one current compound at 105. Each compound 

10 
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represents a cycle ia an iterative process in which multiple compoimds are analyzed to 
generate relevant phenotypic data. Each new iteration begins with operation 105. In 
practice, multiple compounds may be analyzed in parallel, so the iterative/sequential 
nature of the process may not be strictly accurate. Regardless of how the process is 
5 depicted, multiple compounds are evaluated at some point. The flow chart simply 
depicts this fact. 

With a cxarrent compound selected, the process next selects a particular 
combination of compound dose and ceU line for application of the dose. See block 
107. In a preferred embodiment, each compound has an associated matrix of cell lines 

10 and doses. This matrix represents the fact that multiple distinct cell lines are treated 
with the compoimd of interest, each at multiple doses. Each combination of dose and 
ceU liae provides separate phenotypic information. Ultimately, the response curve 
passes through distinct points, each representing a separate dose. At each dose, the 
phenotypic information spans multiple cell lines. In principle, the points on the 

15 response path can be confined to a single cell line. 

After the current combination of dose and cell line has been selected and 
provided, the process next images the cells of the current cell line that have been 
exposed to the current compound at the current dose. See block 109. If more than one 
cell component is to be considered, the imagiag apparatus may generate multiple 

20 images, one for each cell component/maiker combination. At 111, the process 
performs an image analysis that measures and stores parameter values. In some 
embodiments, these features will be separately extracted firom multiple images of the 
cell line taken at different times after exposure to the compound. At 113, the process 
determines whether there are additional combinations of dose and ceU lines to be 

25 considered. If so, process control returns to 107 where the next combination of dose 
and cell line is selected. 

Ultimately, all the relevant combinations of dose and cell line for a given 
compound have been imaged and analyzed. At that point, process control proceeds to 
block 115 where the system combines feature values across multiple cell lines to 
30 obtain separate phenotypic vectors for each separate dose. These phenotypic vectors 
represent the individual points in a response path associated with the current 
compound. 

At 117, the process determines whether there are more compounds to be 

considered as part of the analysis. If so, process control returns to block 105 where 

35 the next current compound is selected. Thereafter, that compoimd is treated as 

described above with respect to blocks 107 through 115. 

11 
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Note that some or all of the operations described above for each compound 
may be automated and performed by a machine. The machine operations may be 
performed by various image acquisition and image analysis apparatus. 

After each of the compounds from the collection identified at 1 03 have been 
5 analyzed as described above, the process has numerous phenotypic vectors 
(quantitative phenotypes), each of which is associated with a particular combination 
of chemical compound and dose. Each of the vectors represents a point in 
multidimensional space. The numerous dimensions may be difficult to depict in a 
manner that presents meaningful information to a human viewer. Therefore, process 

10 101 next reduces the dimensionaUty of feature space and depicts response paths for 
each compound in the reduced dimensional space. See 119. One preferred approach 
to this involves performing a principle component analysis on the collection of 
separate phenotypic vectors. After the reduced dhnensional space has been generated, 
the system may next compare the separate paths of the individual compounds. See 

15 121. This can provide relevant information about the mechanism of action of the 
various compoimds. It allows a human or computer algorithm to compare the various 
paths and draw conclusions about the mechanisms of action of the various 
compoxmds. Note that it is not strictly necessary to depict the response paths in a 
reduced dimensional space prior to comparing the separate paths. Thus, if a 

20 computing device is used to do the comparison, then operation 119 may be optional. 

Note that the discussion of process 101 treats exposure to chemical 
compounds as the stimuli of interest. The process 101 can be extended to cover any 
particular stimulus, not just exposure to chemical compounds. As mentioned, stimuli 
of interest to the present invention include exposure to biological agents, exposure to 
25 various fields, forces, and radiation, deprivation of agents important for normal cell 
growth and fimctioning, etc. 

Also, alternative definitions of response path that do not involve variation over 
dose or time could be employed. For example, a path could be provided through 
multiple distinct cell lines, where each point on tiie path represents a different cell 
30 line. 



Selecting Experiments for ProviH ing Respo nse Paths 

Initially, a relevant collection of stimuli for consideration in the analysis must 
be selected. As mentioned, the stimuH suitable for use with this invention span a wide 
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range of physical agents, forces, fields, etc. Generally, the collection of stinmli 
chosen for a particular analysis may be selected with no prior assumptions. More 
often, the stimuli are selected because they are believed to have related and int^esting 
effects on cells, hi the case of potential therapeutic compoimds, a number of chemical 
5 compoxmds may be selected because they are beUeved to have a similar mechanism of 
action when applied to particular cells. For example, compounds may be selected 
because they are believed to possess anti-mitotic properties when appUed to cancer 
cells. 

The data used for the analysis of the inv^tion may be derived fi'om a wide 
10 range of experiments. Such experiments typically span a matrix of experimental 
conditions. Such matrix may include e^qperimental variations in the choice of 
stimulus, the level of each stimulus, the cell lines to which the stimulus is applied, and 
the particular components within a ceU line that are analyzed. Multiple compoimds 
may be appUed in multiple concentrations to multiple cell lines. For each 
15 combination of compound, dose, and cell line, multiple images may be obtained. 
Each such image contains information about a separate component/marker 
combination within the cell. Note that the invention is not limited to this wide- 
ranging matrix. At its essence, the invention simply involves considering a single 
stimulus at multiple levels. Of course, each such level should provide multivariate 
20 data about a cell phenotype. However, it is unnecessary to employ multiple oell lines 
and/or multiple cellular components in generating the relevant multivariate data. 

The component/marker combinations used in a particular study should be 
chosen based upon the area of interest For example, oncology investigations may 
require a different set of markers than cardiovascular investigations. Further, the 

25 choice of markers should vary over a range of cell biology. For example, it typically 
would be unnecessary to choose two separate markers that both image microtubules. 
DependiQg upon the apphcatidn, the markers can have a very high degree of 
specificity, as in the case of an antibody for tubxilin or can be more lower degree of 
specificity, as in the case of lectins. Note that some lectins, such as Lens culinaris 

30 (LC) lectin actually binds to various polysaccharides. Because most of the time these 
polysaccharides components are enriched in the Golgi, LC lectin still can be an 
effective marker for Golgi. 

Generally, cell components tracked in presently preferable embodiments can 
include proteins, protein modifications, genetically manipulated proteins, exogenous 
35 proteins, enzymatic activities, nucleic acids, lipids, carbohydrates, organic and 
inorganic ion concentrations, sub-cellular structures, organelles, plasma membrane. 
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adhesion complex, ion channels, ion pumps, integral membrane proteins, ceU surface 
receptors, G-protein coupled receptors, tyrosine kinase receptors, nuclear membrane 
receptors, ECM binding complexes, endocytotic machinery, exocytotic machinery, 
lysosomes, peroxisomes, vacuoles, mitochondria, Golgi apparatus, cytoskeletal 
5 filament network, endoplasmic reticulum, nuclear membrane, proteosome apparatus, 
chromatin, nucleolus, cytoplasm, cytoplasmic signaling apparatus, microbe 
specializations and plant specializations. 

The following table illustrates some cell components and markers (labeling 
agents) that may be used in embodiments of the present invention. Other markers can 
10 be used in various embodiments without departing firom the scope of the invention. 



Cell component 


Marker or 


liisease oiaie 




Component 




Plasma membrane 


Carbocyamne dyes 


Apoptosis-Cancer 


(mcluding overall ceU 


Fnospnatiaylsenne 


/\pop to SlS-iN CUlcU 


shape) 


Various lipids 


uegeneranve lys 




Glycoproteins 




Adhesion complexes 


Cadherins 


Thrombosis 




Integrins 


Metastasis 




Occludin 


Woimd healing 




Gap junction 


Xnflammatory Ds 




ERM proteins 


Dermatologic Ds 




CAMS 






Catenins 






Desmosomes 




Ion Channels and Pumps 


Na/K Atpase 


Cystic fibrosis 




Calcimn channels 


Depression 




Serotonin reuptake pump 


Congestive Heart Failure 




CFTR 


Epilepsy 




SERCA 




G coupled receptors 


P adrenergic receptor 


Hypertension 




Angiotensin receptor 


Heart Failiure 






Angina 
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Tyrosine kinase receptors 


PDGF receptor 
FGF receptor 
IGF receptor 


Cancer 

Wound healing 
Angiogenesis 
Cerebrovascular Ds 


ECM binding complexes 


Dystroglycan 
Syndecan 


Muscular Dystrophy 


Endocytotic machinery 

1 


Clathrin 

Adaptor proteins 
COPs 
Presenilins 
Dynamin 


Alzheimer*s Ds 


Exocytotic machinery 


SNAREs 
Vesicles 


Epilepsy 
Tetanus 

SyRtemin TnflaTmnatioTi 

Allergic Reactions 


Lysosomes 


Acid phosphatase 
Transferrin 
Lysotracker Red 


Viral diseases 


PeroxisomesA/^ acuoles 




Neural degenerative Ds 


Mitochondria 


Caspases 

Apoptosis inducing factor 
Fl ATPase 
Fluorescein 
Cyclo-oxygenase 
Mitotracker Red 
Mitotracker Green 


Apoptosis 

Neural degeanerative Ds 
Mitochondrial Cytopathies 
Inflammatory Ds 
MetaboHcDs 


Golgi Apparatus 


Lens culinaris lectin 
DiOC6 carbocyanine dye 
COPs 

Antibodies specific for Golgi 
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Cytoskeletal 
Networks 



Filament j Microtubxiles 
Actin 

Intermediate Filaments 
Kinesin, dynein, myosin 
Microtubule associated 
proteios 

Actin binding proteins 
Rac/Rho 
Keratins 
GFAP 

Von Wiltbrand's factor 



Cancer 

Neural degenerative Ds 
Inflammatory Ds 
Cardiovascular Ds 
SldnDs 



Endoplasmic Reticidmn 



SNARE 
PDI 

Ribosomes 



Neural degenerative Ds 



Nuclear Membrane 



Lamins 

Nuclear Pore Complex 



Cancer 



Proteosome Apparatus 
Chromatin 



Ubiquityl transferases 



Cancer 



Nucleolus 



DNA 

Histone proteins 
Histone deacetylases 

Telomerases 

Phase markers 



Cancer 
Aging 



Cytoplasm 



Intermediary 

Enzymes 

BRCAl 



Metabolic Cancer 



Cytoplasmic 
Apparatus 



Signaling | Calcium 
Camp 
PKC 





Cardiovascular Ds 
Migraine 
Apoptosis 
Cancer 



Microbe Specializations 



Flagella 
Cilia 

Cell Wall components: 
Chitin synthase 



Infectious Ds 



Plant specializations 



Choloroplast 

Cell Wall components 



Crop Protection 
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In one preferred embodiment, the cellular components considered in separate 
images include one ore more of DNA, cytoskeletal proteins, and Golgi. In a specific 
embodiment, the images for each combination of cell line, dose, and compound 
5 include a DNA image, a tubulin image, and a Golgi image. Various markers can be 
used for each of these components. In a preferred embodiment, the DNA marker is 
DAPI, the tubulin marker is an antibody specific for tubulin, and the Golgi marker is 
LC lectin. 

In one specific approach, the above three markers are analyzed using two 
10 separate processes. In a first process, a cell line is simply stained with a marker for 
DNA. In a second run, the cell line is stained with all three markers. The first process 
run is used to simply identify cell cycle information. For example, this run is used to 
determine the proportion of cells in each separate phase of the cell cycle (Gl, S, G2, 
M, and/or various subphases of M). The two process runs are employed because 
15 imaging tubulin and Golgi require repeating washing of the cells. This process 
selectively causes some cells to wash away; specifically rounded up and mitotic cells. 
Therefore, the remaining cells imaged for tubuhn and Golgi are biased toward 
interphase states. 

Regarding the doses or "levels" of the various stimuli, one should endeavor to 
20 choose a range of doses that define and active zone for affecting phenotype in a cell 
line of interest. In one approach, researchers perform a preliminary experiment with 
each drug. The preliminary experiment may involve titration across a wide range of 
concentrations. The titration may measure cell count or other appropriate biological 
parameter. An upper boundary of the active zone may be a concentration at which 
25 further increases of concentration have no additional affect on the cells. For example, 
the upper boundary may be the minimum concentration at which all cells are killed. A 
lower bound of the active zone is the lowest concentration at which some biological 
affect can be observed. 

In some cases, the highest dose allowed by the process is governed by some 
30 physical parameter such as the maximum solubihty of a compoimd. Alternatively, it 
may be governed by the maximum volume of a compoimd solution that can be 
administered to a well without having the solvent significantly affect the cells. 

In a preferred embodiment, a highest level of the stimuli is first identified by 
some technique. Then, additional lower level of the stunulus are identified by 
35 incremental reductions. For example, in the case of a chemical compoimd, serial 
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dilutions may be performed to generate lower level doses. At a miiumiim, at least two 
levels of the stimulus must be considered. Preferably significantly more levels are 
considered. In a preferred embodiment, at least five separate stimulus levels are 
considered. In a specific preferred embodiment, eight separate levels are considered. 
5 If a chemical compoimd serves as the stimulus, then the highest concentration of the 
compound should be at least about two times that of the lowest concentration. 

As indicated, phenotypic vectors for given stimuU and level combinations may 
include multivariate information taken fi*om different cell lines. However, this need 
not be the case, as all the muWvariate data of interest may be obtained &om a single 

10 cell line. Generally, a researcher will chose one or a range of cell lines that are 
relevant to the area of interest. For example, if the researcher focuses on oncology 
applications, the cell lines chosen may include different types of cancers and possibly 
other cells lines that allow one to identify typical side effects of anti cancer dmgs. In 
one specific embodiment pertaining to oncology, six different cell lines are 

15 considered. These include HUVEC (human umbilical vein endothelial cells), A498, 
A548, SF268, SKOV3, andDU145. 



Imaging 

As indicated, the phenotypic data characterizing each point on a response 
20 curve is derived, at least in part, from images of cell lines exposed to particular 
combinations of stimulus type and stimulus level. See block 109 in Figure 1, for 
example. Various techniques for preparing and imaging appropriately treated cells are 
described in U.S. Patent Applications 09/310,879, 09/311,996, and 09/311,890, 
previously incorporated by reference. In the case of cells treated with a fluorescent 
25 marker, a collection of such cells is illuminated with Ught at an excitation frequency. 
A detector is tuned to collect light at an emission frequency. The collected light is 
used to generate an image, which highlights regions of high marker concentration. 

Additional operations may be performed prior to, during, or after the imaging 
operation (109) of Figure 1. For example, "quahty control algorithms" may be 
30 employed to discard image data based on, for example, poor exposure, focus failwes, 
foreign objects, and other imaging failures. Generally, problem images can be 
identified by abnormal intensities and/or spatial statistics. 

In a specific embodiment, a correction algorithm may be appUed prior to 
segmentation to correct for changing light conditions, positions of wells, etc. In one 
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example, a noise reduction technique such as median filtering is employed. Then a 
correction for spatial diflEerences in intensity may be employed. In one example, the 
spatial correction comprises a separate model for each image (or group of images). 
These models may be generated by separately summing or averaging all pixel values 
5 in the x-direction for each value of y and then separately summing or averaging aH 
pixel values in the y direction for each value of x. In this manner, a paraboUc set of 
correction values is generated for the image or images \mder consideration. Applying 

■ 

the correction values to the image adjusts for optical system non-linearities, mis- 
positioning of wells during imaging, etc. 

10 The production of the images includes cell plating, compound dilution, 

compoxmd addition and imaging focusing. Failures in any these systems can be 
detected by a variety of methods. For example, cell plating could fail because of a 
clogged tip in a delivery pipette. Such failure can be identified by adding a 
fluorescent dye or bead to the cell suspension. The fluorescence of this dye or bead is 

15 chosen to be at a different channel (wavelength) than the markers used to image 
cellular components. Another potential failure could occur during compound 
delivery. To detect such failures, one can add a fluorescent dye or bead in the 
compoimd plate before compoimd dilution. The amoxmt of fluorescent dye or bead is 
proportional to the amoimt of compoxmd. Yet another potential problem occurs when 

20 the focus of the image acquisition system changes during imaging. To account for 
such spatial biases, one can employ control wells containing, for example, cells with 
no or neutral compounds interspersed throughout the plate. Still another problem 
results from foreign objects (e.g., small dust particles) in the well. This can be 
addressed with image segmentation and statistical outlier identification techniques. 

25 Generally the images xised as the starting point for the methods of this 

invention are obtained from cells that have been specially treated and/or imaged under 
conditions that contrast the cell's marked components from other cellular components 
and the background of the image. Typically, the cells are fixed and then treated with a 
material that binds to the components of interest and shows up in an image (i.e., the 

30 marker). Preferably, the chosen agent specifically binds to DNA, but not to most 
other cellular biomolecules. 

Multivariate Phenotvpic Data from Images 

At every combination of dose, cell line, and compound, one or more images 
35 are obtained. As mentioned, these images are used to extract various parameter values 
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of relevance to a biological, phenotypic characterization of the compound of interest. 
Generally a given image of a cell, as represented by one or more markers, can be 
analyzed to obtain any number of image parameters. These parameters are typically 
statistical or morphological in nature. The statistical parameters typically pertain to a 
5 concentration or intensity distribution or histogram. The parameters chosen for use 
with this invention should relate to the expected biological response of the cell lines to 
the compound of interest. The parameters should also represent a diverse set of 
phenotypic characteristics. 

Some general parameter types suitable for use with this invention include a 
10 cell count, an area, a perimeter, a length, a breadth, a fiber length, a fiber breadth, a 
shape factor, a elliptical form factor, an inner radius, an outer radius, a mean radius, 
an equivalent radixxs, an eqitivalent sphere volume, an equivalent prolate volume, an 
equivalent oblate volimie, an equivalent sphere surface area, an average intensity, a 
total intensity, an optical density, a radial dispersion, and a texture diflference. These 
15 parameters can be average or standard deviation values, or frequency statistics fix>m 
the descriptors collected across a population of cells. In some embodiments, the 
parameters include features from different cell portions or cell types. 

Examples of some specific parameters/descriptors that may be suitable for use 
in multivariate response paths of this invention are included in the following table. 
20 Other descriptors can also be used without departing from the scope of the invention. 



Name of Parameter 


Explanation/Comments 


Count 


Number of objects 


Area 




Perimeter 




Length 


X axis 


Widfli 


Y axis 


Shape Factor 


Measure of roundness of an object 


Height 


Z axis 


Radius 




Distribution of Brightness 




Radius of Dispersion 


Measure of how dispersed the marker is from its 
centroid 


Centroid location 


x-y position of center of mass 


Nimiber of holes in closed objects 


Derivatives of this measurement might include, for 
example, Euler number (— number of objects - 
nimiber of holes) 


Elliptical Fourier Analysis (EFA) 


Multiple frequencies that describe the shape of a 
closed object 


Wavelet Analysis 


As in EFA, but using wavelet transform 
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Interobject Orientation 


Polar Coordinate analysis of relative location 


Distribution Interobiect Distances 


Including statistical characteristics 


Spectral Output 


Measures the wavelength spectrum of the reporter 
dye. Includes FRET 


Optical density 


Absorbance of light 


Phase density 


Phase shifting of light 


Reflection interference 


Measure of the distance of the cell membrane from 
the surface of the substrate 


1 9 and dimensional Fourier 


Spatial frequency analysis of non closed objects 


1,2 and 3 dimensional Wavelet 

Arifilv<5iQ 


Spatial frequency analysis of non closed objects 


Eccentricity 


The eccentricity of the ellipse that has the same 
second moments as the reeion. 
A measure of object elongation. 


Long axis/Short Axis Length 


Another measure of object elongation. 




Perimeter of the smallest convex polygon 
surrounding an object 


r^nnvex area 


Area of the smallest convex polygon surrounding an 
obiect 


Solidity 


Ratio of polygon bounding box area to object area. 




proportion of pixels in the bounding box that are also 
in the region 


Granularity 




Pattern matching 


Significance of similarity to reference pattern 


Volume measurements 


As above, but adding a z axis 


Number of Nodes 


The number of nodes protruding from a closed 
obiect such as a cell; characterizes cell shape 


End Points 


Relative positions of nodes from above 



The features used in the actual points comprising a response path of this 
invention may be parameters directly extracted from the images or they may be 
biological characterizations derived from the parameters. Note that the points may 
also mclude some features that were not directly or indirectly obtained from the 
5 images. For example, the points may include information obtained from public 
sources such as databases, literature, etc. Further, the features comprising the points 
may include non-image related data, such as data obtained from chemical and 
biological assays. 

Often, the parameters are chosen based upon a biological imderstanding. For 
10 example, if a cell's state in the cell cycle is important to the biological problem being 
investigated, then parameters that characterize the amoamt of DNA in a cell and/or the 
degree of condensation of that DNA into chromosomes is relevant. In a specific 
example, cell cycle parameters include the total quantity of DNA in a nucleus, the area 
of the nucleus, and the intensity variance of the cellular DNA. A fiill discussion of the 
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relevant parameters for characterizing the cell cycle is presented in U.S. Patent 
Application No.09/729,754, previously incorporated by reference. 

Similarly, if an objective is to characterize the Golgi in a cell, this can be 
accomplished with parameters that define the location of the Golgi with respect to the 
5 nucleus, describe the texture of the Golgi and describe the local concentration of 
Golgi components. The fiill discussion of the parameters relevant in characterizing 
Golgi is presented m U.S. Patent AppUcation No. 09/792,012 (Atty. Docket No. 
CYTOP012), previously incorporated by reference. Of specific interest, the Golgi 
complex in a perinuclear region may be characterized using parameters such as the 
10 mean, standard deviation, and kurtosis of pixel intensity, and various eigenvalues 
obtained by singular value decomposition of a pixel intensity matrix for the Golgi 
marker. From these parameters, the Golgi complex of a given cell may be 
characterized as normal, difiused, dispersed, or dispersed and diffused. 

Further, if the cell shape provides further relevant phenotypic data, then 
15 parameters can be chosen accordingly. In one embodiment, a tubuhn or other 
cytoskeletal component is marked and imaged to provide parameters relevant to cell 
shape. Specific examples of such parameters include the number of nodes on a cell 
image, the distance between end point of those nodes, a coefficient of tubulin 
polymerization (e.g., average pixel intensity of object pixels in a tubulin channel), 
20 averaged across all cells in a population, and a coefficient of microtubule 
reorganization (e.g., standard deviation of wavelet coefficients), averaged across all 
cells in a population. A fiiU discussion of parameters relevant to characterizing ceU 
shape can be found in U.S. Patent Application No. 09/792,013 (Atty. Docket No. 
CYTOPO 13), previously incorporated by reference. 

25 While a fiindamental biological understanding can often direct one to the 

appropriate choice of parameters for use ia this invention, a systematic analysis of 
data can help identify parameters that might not be immediately apparent. Such 
analysis can be conducted in a mamier that finds parameters that are best able to show 
subtle differences in the response path. By considering the effect of varying a single 

30 parameter at a time, one can quickly hpme in the most relevant parameters for 
developing response curves in accordance v/ith this invention. 

As applied to specific markers, one preferred collection of parameters includes 
(1) total nxnnber of cells in a population of interest, (2) niraiber of cells relative to 
nimiber of cells in one or more controls, (3) proportion of cells in each of the stages of 
35 cell cycle (Gl, S, G2, pre-anaphase mitotic, post-anaphase mitotic), (4) area of cell 
nuclei, averaged for each of the stages of the cell cycle, (5) diameter of cell nuclei, 
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averaged for each of the stages of cell cycle, (6) axes ratio (measure of elongation) of 
nuclei, averaged for each of the stages of cell cycle, (7) eccentricity of nuclei (another 
measure of elongation), averaged for each of the stages of cell cycle, (8) solidity of 
nuclei (measure of a shape), averaged for each of the stages of cell cycle, (9) total 
5 intensity of nuclei pixels (a measure of amount of DNA in a cell), averaged for each 
of the stages of cell cycle, (10) variance of pixel intensities in a nucleus, averaged for 
each of the stages of cell cycle, (11) proportion of cells with normal Golgi in a 
subpopulation of cells in each stage of the cell cycle, (12) proportion of cells with 
diffuse Golgi in a subpopulation of cells in each stage of the cell cycle, (13) 

10 proportion of cells with dispersed Golgi in a subpopulation of cells in each stage of 
the cell cycle, (14) proportion of cells with dispersed aud diffused Golgi in a 
subpopulation of cells in each stage of the cell cycle, (15) coefficient of the Golgi 
dispersion (kurtosis of the Golgi marker in a region of the cell), (16) coefficient of 
tubulin polymerization (average pixel intensity of object pixels in a tubulin chaimel), 

15 averaged across all cells in a population, and (17) coeflBcient of microtubide 
reorganization (standard deviation of wavelet coefficients), averaged across all cells in 
a population. 

In a more specific embodiment, the parameter set for each dose/compo\md 
point in a response path includes parameters derived fi*om DNA markers, Golgi 
20 markers, and tubuUn markers. It is been found that the foUowdng ten parameters 
provide particularly useful phenotypic results. 

1 . The size of the nuclei as derived firom a DNA marker. This value is 
provided as the average area of the nuclei of aU interphase cells in an image. 

2. The average elUpsicity of nuclei of interphase cells. 

25 3. The difference between the proportion of interphase cells in Gi and 

the proportion of interphase cells in G2. 

4. The proportion of interphase cells in S. 

5. The "mitotic index" which specifies the proportion of mitotic cells 

in the image. 

30 6, The proportion of interphase cells having normal Golgi. 

7. The proportion of interphase cells having diffuse Golgi. 
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8. The proportion of interphase cells having diffuse and disperse 

Golgi. 

9. The mean pixel intensity obtained from marked tubulin. This value 
is obtained by determining an overall threshold for the image and considering only 

5 objects having intensities above that threshold. The mean value of intrasity of aU 
pixels over the threshold is used. 

10. A second order wavelet obtained from all pixels over the threshold 
value in the tubulin image. 

In a particularly preferred specific embodiment, the above ten parameters are 
10 obtained for each of multiple cell lines, such as the six cell lines defined above, and 
packaged in a vector. This vector then represents a single dose/compound point to be 
used in a response path for the compound of interest. Note that this vector spans 
multiple ceU lines and multiple phenotypic features. To obtain these various features, 
three separate markers were considCTed. A DNA marker was used for the first five 
15 ^ parameters. A Golgi marker was used for the next three parameters. Finally, a tubulin 
marker was used for the last two parameters. 

VisuaUzation and Comparison of Response Paths 

Lists of stimuh and associated quantitative phenotypes may be stored as 
20 database records or other data structures that can be queried or otherwise accessed as 
part of an analysis procedure. The stimuh may also be associated with other relevant 
data such as clinical toxicity, cellular toxicity, hypersensitivity, mechanism of action, 
etc. (when available). The stored phenotypic data is used to generate and depict 
response paths. 

25 Various techniques may be employed to visualize the response paths generated 

as described above. In order for a human observer to make meaningful comparisons, 
the space in which the response paths are presented should be comprehendible. Note 
that for compUcated quantitative phenotypes representing individual points on the 
path, there may be very many separate variables (60 in the example above). In 

30 principle, each of these variables represents a separate dimension. So one may be 
confronted with a 60 dimensional space, for example. Obviously, it becomes difficult 
to visualize meardngfiil trends or clusters in high dimensional space. Consider the 
problem of trying to visualize a trend in phenotypes comprised of three cellular 
components (e.g. tubulin, DNA and Golgi), each of which has multiple relevant 
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parameters (e.g. total quantity of DNA and variance in the concentration of DNA). 
Obviously, there are more than three relevant dimensions to be considered in 
analjCTng such phenotypes. 

One possible solution to the problem involves selecting two or three 
5 dimensions (features) that are expected to be most relevant to a particular response 
curve. Unfortunately, it becomes impossible to view other potentially meaningful 
phenotypic features on the same two or three-^dimensional space. 

Various techniques may be employed to address this problem. Such 
techniques create a lower dimensional space in which the individual dimensions 

10 capture two or more features of the data. Examples of such techniques include 
principle component analysis, linear and non-linear discriminant analysis, 
multidimensional scaling, and projection pursuit techniques. A particularly preferred 
approach involves the use of principle component analysis. Principle component 
analysis determines the vectors (dimensions) through which a data set shows the 

15 greatest variation in multidimensional space. The first principle component shows the 
direction of greatest variation in the data. The second principle component shows the 
direction of the second greatest variation in data and so on. One can select as many 
principle components as are suitable to depict one's data. Typically, the first one, two, 
or three principle components are selected for presenting data to human observers. 

20 Principal component analysis is described more fully in Jackson, J. E. (1991) A User 
Guide to Principal Components. New York: John Wiley and Sons; and Jolliffe, I. T. 
(1986) Principal Component Analysis. New York: Springer- Verlag, both of which are 
incorporated herein by reference for all purposes. 

Various commercially available tools for performing principle component 
25 analysis are available. One suitable statistical computiag package for performing 
PGA is available fi-om Insightfiil Corporation (formerly MathSofl) of Seattle, WA. 
Principal component analysis can be applied to quantitative phenotypic data sets in a 
straight-forward manner. However, it will generally be necessary to standardize 
phenotypic data sets before submitting them to principle component analj^is. This is 
30 because the various scalars that comprise the individual features of a quantitative 
phenotype reside on vastly different scales. For example, the mitotic index will range 
firom zero to one hxmdred percent, while the size of the nuclei, average ellipsicity of 
the nuclei, average pixel intensity of the tubulin marker, etc. each have very different 
scales and associated imits. To bring these various featm*es onto a comparable scale 
35 for meaningful PCA analysis, one may perform transformations to standardize the 
data. In one preferred embodiment, each of the dimensions is scaled by considering 



25 



wo 02/067182 



PCTAJS02/05553 



all the data along that dimension (e.g., all values of nucleus area), subtracting the 
mean of that data and dividing by the standard deviation. This effectively scales the 
data for standardization. 

After PCA is performed by a suitable tool, the results should be presented 
5 graphically. Various graphical tools are suitable for this purpose. One is provided by 
S+ Corporation. Another particularly useful graphical depiction tool is Spotfire.net 
available firom SpotFire, hic. of Cambridge, MA. Any of these tools will not only 
present the data in principle component space, but will also identify which variables 
(features) contribute the most to each of the principle components. Because the 
10 present invention is concerned primarily with paths, the graphical depiction will 
preferably show connections between individual points along the dose response path 
for each particular stimulus. 

Meaningful comparison between related stimuli requires that one identify 
patterns or trends in various response curves. This may be accompUshed with or 

15 without the aid of a visualization technique/tool of the type described above. If a 
human observer is to participate in the pattern recognition, then such a visualization 
tool will typically provide great assistance. However, if a machine is to do the 
comparison, the reduced dimensional visualization technique may be unnecessary. 
Examples of techniques that may be employed for such comparison include 

20 techniques that determine an average difference or distance between two potentially 
related stimulus response curves, clustering, and the like. 

Various examples applying principle component analysis to complex 
multivariate quantitative phenotypes are presented below. In these examples, dose 
response curves through multivariate phenotypes clearly show trends and clustering 
25 based upon mechanism of action. 

In comparing response paths of potentially related stimiili, one typically 
identifies similarities in the general pathways. These similarities often show up in the 
path trajectories, starting and ending points, etc. In comparing stimuli such as 
exposure to chemical compounds, note that most drugs have response paths that 
30 follow reasonably similar trajectories at relatively low concentrations, but then diverge 
as the concentration increases. In other words, the phenotypic manifestations of 
distinct mechanisms of action only appear most pronounced at high concentrations. 
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Software/Hardware 

Generally, embodiinCT.ts of the present invention employ various processes 
involving data stored in or transferred through one or more computer systems. 
Embodiments of the present invention also relate to an apparatus for performing these 
5 operations. This apparatus may be specially constructed for the required purposes, or 
it may be a general-purpose computer selectively activated or reconfigured by a 
computer program and/or data structure stored in the computer. The processes 
presented herein are not inherently related to any particular computer or other 
apparatus. In particular, varioxis general-purpose machines may be tised with 
10 programs written in accordance with the teachings herein, or it may be more 
convenient to construct a more specialized apparatus to perform the required method 
steps. A particular structure for a variety of these machines will appear from the 
description given below. 

In addition, embodiments of the present iaveation relate to computer readable 
15 media or computer program products that include program instructions and/or data 
(including data stmctures) for performing various computer-implemented operations. 
Examples of computer-readable media include, but are not hmited to, magnetic media 
such as hard disks, floppy disks, and magnetic tape; optical media such as CD-ROM 
disks; magneto-optical media; semiconductor memory devices, and hardware devices 
20 that are specially configured to store and perform program instructions, such as read- 
only memory devices (ROM) and random access memory (RAM). The data and 
program instmctions of this invention may also be embodied on a carrier wave or 
other transport medium. Examples of program instructions include both machine 
code, such as produced by a compiler, and files containing higher level code that may 
25 be executed by the computer using an interpreter. 

Figure 2 illustrates a typical computer system that, when appropriately 
configured or designed, can serve as an image analysis apparatus of this invention. 
The computer system 200 includes any nvmiber of processors 202 (also referred to as 
central processing units, or CPUs) that are coupled to storage devices including 

30 primary storage 206 (typically a random access memory, or RAM), primary storage 
204 (typically a read only memory, or ROM). CPU 202 may be of various types 
including microcontrollers and microprocessors such as programmable devices (e.g., 
CPLDs and FPGAs) and unprogrammable devices such as gate array ASICs or general 
piirpose microprocessors. As is well known in the art, primary storage 204 acts to 

35 transfer data and instructions uni-directionally to the CPU and primary storage 206 is 
used typically to transfer data and instructions in a bi-directional manner. Both of 
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these primary storage devices may include any suitable computer-readable media such 
as those described above. A mass storage device 208 is also coupled bi-directionally 
to CPU 202 and provides additional data storage capacity and may include any of the 
computer-readable media described above. Mass storage device 208 may be used to 
5 store programs, data and the like and is typically a secondary storage medium such as 
a hard disk. It will be appreciated that the information retained within the mass 
storage device 208, may, in appropriate cases, be incorporated in standard fashion as 
part of primary storage 206 as virtual memory. A specific mass storage device such as 
a CD-ROM 2 14 may also pass data uni-directionally to the CPU. 

10 CPU 202 is also coupled to an interface 210 that connects to one or more 

input/output devices such as such as video monitors, track balls, mice, keyboards, 
microphones, touch-sensitive displays, transducer card readers, magnetic or paper tape 
readers, tablets, styluses, voice or handwriting recognizers, or other well-known input 
devices such as, of course, other computers. Finally, CPU 202 optionally may be 

15 coupled to an extemal device such as a database or a computer or telecommunications 
network using an extemal connection as shown generally at 212. With such a 
connection, it is contemplated that the CPU might receive information from the 
network, or might output information to the network in the course of performing the 
method steps described herein. 

20 In one embodiment, the computer system 200 is directly coupled to an image 

acquisition system such as an optical imaging system that captures images of cells. 
Digital images from the image generating system are provided via interface 212 for 
image analysis by system 200. Alternatively, the images processed by system 200 are 
provided from an image storage source such as a database or other repository of cell 

25 images. Again, the images are provided via interface 212. Once in the image analysis 
apparatus 200, a memory device such as primary storage 206 or mass storage 208 
buffers or stores, at least temporarily, digital images of the cells. In addition, the 
memory device may store the quantitative phenotypes that represent the points on the 
response path. The memory may also store various routines and/or programs for 

30 analyzing the presenting the data, including the response paths. Such 
programs/routines may include programs for performing principal component 
analysis, regression analyses, path comparisons, and for graphically presenting the 
response paths. 
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Examples 

As indicated above, an underlying premise of this invention is that changes in 
cell physiology are detectable flirough single-cell image analysis of cellular markers 
and, by measuring such changes in a sophisticated fashion, one can monitor biological 
5 similarities and differences between compounds. To characterize and demonstrate the 
ability of multivariate quantitative phenotype response paths to distinguish different 
classes of molecules, a panel of commercially available compounds with broad 
mechanisms of action was tested using the present invention. The types of 
compotmds tested are listed in the following table. 



Class 


Primary Protein Target 


- 

Number of 
Compouods 


Calcium 


Endoplasmic reticulum Ca2+- 
ATPase 


3 


Calcium 


Calmodulin 


6 


Cytoskeleton 


Actin 


8 


Cytoskeleton 


Tubulin 


9 


G Protein Effectors 


G-proteins G,- and Go 


4 


Gene Regulation 


Topoisomerase n 


6 


Ion Pump 


V-ATPase 


2 


Oxidative Phosphorylation 


Mitochondrial ATPases 


3 


Posttranslational 
Modification 


Famesyltransferase 


2 


Posttranslational 
Modification 


Geranylgeranyltransferase 1 


3 


Protein Kinase 


p38 MAP kinase 


3 


Protein Kinase 


PKC 


3 


Protein Kinase 


p34cdc2/cyclin B 


4 



A three-dimensional representation of three principal components calculated 
for each compound is shown in Figure 3A. The data points are presented within the 
first three principal components obtained for the entire data set. Note that the controls 

1 5 cluster tightly in the center of the graph. The specific features and cell lines used to 
constract the quantitative phenotypes shown as points ia the plot are set forth above in 
the '^Multivariate Phenotypic Data from Images" and "Selecting Experiments for 
Providing Response Paths" sections. Specifically, the features are the ten features 
listed near the end of the "Multivariate Phenotypic Data from Images" section 

20 collected for each of the six cell lines hsted in the "Selecting Experiments for 
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Providing Response Paths" section. In Figure 3A, the connected data points link 
identical compounds at increasing concentration. Specifically, there were eight 
different concentrations for each compound. 

This dose-response path has proven to accurately classify different compound 
5 mechanisms. Even though maay of the compound classes tested do not have a direct 
effect on the components/markers used in this experiment, quantitative phenotype 
paths are able to detect differences in this wide array of compounds. For example, 
there is a distinct separation of the famesyltransferase and geranylgeranyltransferase 
inhibitors from each other even though both are involved in post-translational protein 
10 modificatioiL These compounds are quite different from compounds that inhibit 
mitochondrial function (see Figure 3B). 

Similarly, actin and tubulin inhibitors show distinct dose-response trajectories 
overall even though both classes affect the cytoskeleton (see Figure 3C). 
Furthermore, compoimds that effect signaling pathways directly, such as protein 
15 kinase inhibitors and calcium sensitizes, are uniquely differentiated from the other 
classes as well (see Figure 3D). 

The demonstration that many compound classes follow unique dose-response 
trajectories validates this invention's ability to measure the inhibition of many 
different types of targets by classifying a quantitative phenotype directly using a few 
20 carefully selected biological analyses. This information may be used to prioritize and 
expand marker sets, cell lines, and time-points used to generate the quantitative 
phenotypes. Further, it greatly increases the abihty to resolve and differentiate ever 
more similar compound mechanisms. 

Another example demonstrates the ability to investigate biological feature data 
25 (of the type used in the first example) to understand the differences and similarities 
between compoimds. Despite the fact that the marker sets in this example did not 
include components of the actin cytoskeleton, the actin inhibitors are all uniquely 
classified by the quantitative phenotype paths. Figure 4 highlights the deviation of the 
actin inhibitor Cytochalasin A from the other actin inhibitors. Figure 4 also shows 
30 that the Cytochalasin A deviation is reliably reproduced when this experiment was run 
a second time. 

To identify which biological features were changing in Cytochalasin A as the 
concentration of compound increased, biological feature image plots (not shown) for 
all cell lines and concentrations of two actin inhibitors were compared. The most 
35 significant difference in these image plots is the lack of an increase in the tubulin 

30 
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features at increasing concentrations of Cytochalasin A. All of the other actin 
inhibitors showed an increase in the tubulin biological features at the same 
concentrations (data not shown). 

The change in the tubulin feature was fiirther validated by inspection of the 
5 dose response graphs of all actin inhibitors for that feature across all cell lines and 
concentrations (see Figure 5 A (Cytochalasin A) and Figure 5B (Cytochalasin J)) and 
inspection of a representative dose response image montages (not shown). (N.B. 18 
different montages are available for 6 cell lines by 3 replicates.) Review of the 
Uterature confiraied that Cytochalasin A iaterferes with microtubule assembly by 
10 reacting with sulfhydryl groups. The data obtained with the present invention 
suggests that this side effect has a lower affinity than Cytochalasin A has for actin 
itself, as evident from the common PCA trajectory it shares with the rest of the actin 
iohibitors at lower concentrations. 

Another way to review the data is to ask how similar the dose response paths 
15 are for different cell lines. These cell lines, with their different expression patterns, 
may exhibit increased or decreased sensitivity or off-target effects to different 
members within a single class of compoxmds. Because PCA analysis is only shedding 
light on the features that most distinguish compounds, an in depth cell-line sensitivity 
analysis may reveal more subtle differences and insight into molecular specificity. 

20 A way to visuaUze the similarities of compoimds at biologically relevant 

concentration is to vise hierarchical clustering. A dendrogram , in Figure 6 shows the 
compounds profiled at the IC50 for A549 cells. This visualization presents a different 
insight into compoimd mechanism than the PCA plot. For one, the p38cdc2/cyclinB 
inhibitors, the tubulin depolymerizers, and the G-protein activators are all highly 

25 correlated, while differences between the actin inhibitor molecules are highlighted. 
The biological features that are different can be visualized using an image plot of the 
same compounds in cluster order. This type of comparison can shed even more light 
into the biological similarities and differences between compounds at an even more 
subtle level. 

30 The effects of cell line sensitivity are apparent when one inspects the dose 

response graphs for a biological feature for all cell lines (see Figure 7A). Mastoparan, 
an amphiphic wasp poison known to activate G proteins, and its analog 
tetradecapeptide derivative, MAS 7 have very similar quantitative phenotype dose 
response profiles, as seen in the PCA plots (see Figure 3D) and dendrogram (see 

35 Figure 6) yet quite distinct cell line sensitivities (Figure 7B). MAS 7 is reported to 
have 5-fold greater potency than Mastoparan, but the dose response curves show that 
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MAS7 is a more potent compoimd in flie cancer cell lines; it does not exhibit 
increased sensitivity in the normal cell line. Further, the increase in potency is cell 
line dependent. The graphs of Figure 7 A and 7B depict the decrease in the numb&r of 
cells per image for each cell line as a function of concentration. 

5 The oncology program Cytokinetics, Inc., South San Francisco, CA used the 

quantitative phenotyping technologies of this invention in a retrospective study to 
measure cellular phenotypes and biological effects of compounds that inhibit 
oncology program targets. In this effort, Cytokinetics' scientists submitted a series of 
57 primary hits and optimized compounds against Oncology targets for profiling with 
10 this invention. The quantitative phenotyping of this invention identified a numb^ of 
attractive chemical series against a validated target with a defined cellular 
morphological change expected from inhibiting the target. Quantitative phenotyping 
also identified compounds that exhibited off-target biological effects. 

Quantitative profiles of the oncology compounds were compared to a number 
15 of control compoxmds. Inspection of the quantitative profiles in principal conoponent 
space revealed a significant and reproducible separation among the oncology targets, 
non-specific compounds, and controls. Notably, quantitative phenotyping 
differentiated compounds with similar biological effects, but different targets. Figure 
8 shows the dose response path (obtained as described in the examples above) taken 
20 by 83 different compounds clustered in groups representative of the mechanism of 
action. These paths connect identical compounds at increasing concentrations. Note 
that compounds that inhibit the Cytokinetics target are distinct from those that affect 
other cancer targets such as topoisomerase II, tubulin stabilizers, or tubulin 
depolymerizers. 

25 Most of the primary hits and early analogs from multiple chemical classes for 

one of the targets line up along a similar dose-response trajectory, while a few veer 
away (see Figures 9 and 10). This level of multivariate data group structurally 
unrelated compoxmds into groups with similar cellular phenotypes. Having this data 
available early in the drug discovery process allows the research scientist to know if 

30 stracturally imrelated compounds cause a similar cellular phenotype which correlate to 
similar target selectivity, and to rapidly identify and reject compoimds which correlate 
to off-target effects. The more compounds screened with the technologies of this 
invention, the clearer the distinctions become. 

The optimized oncology compoimds superimposed on the dose-response 
35 trajectory of the primary hits and formed a compact cluster of data points in one 
region of the graph that excluded other compoimds (see Figure 11). Inspection of the 
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data and images revealed that the compounds in this region caused the morphological 
phenotype expected from inhibiting the target. 

This region was magnified to display each compound and concentration that 
caused the optimal profile (see Figure 12). There were 39 compounds at many 
5 different concentrations showing that many of the chemical analogs have a similar 
profile. All compounds in the cluster are at concentrations lower than 10 micromolar. 
As one of the criteria for lead componnd selection was potency, the number of 
compoimds was fiu1;her reduced to 8 by asking which compounds were able to 
produce the target specific profile at a delivered concentration of less than 40 
10 nanomolar. 

With the reduced data set, the biological feature data was further inspected to 
gain insight into the similarities and differences between these compounds. The 
compounds selected had a very distinct set of biological features that would have been 
almost impossible to predict in any other way. The differences between these 
15 compounds were fiarther evaluated on a feature-by-feature level to provide fijrther 
insight into the compoimds' mechanism. 

Figure 13 shows how a single multivariate dose response experiment was used 
to identify and narrow down 57 hits to eight higjily potent and specific compounds. 
This demonstrates the breadth of mformation generated by present invention. 

20 As a postscript, of the eight compounds shown in this example, four had been 

selected for extensive toxicity testing, and a derivative of one of those four was 
selected as a development compound. Going forward, the oncology program can use 
the profile information to select backup candidates for this target either by exploring 
the other four compounds identified, or profiling new compoimds, and asking how 

25 similar they are to the optimal profile identified in this experiment. 



Conclusion 

Although the above has generally described the present invention according to 
specific processes and apparatus, the present invention has a much broader range of 
30 applicabiUty. In particular, the present invention has been described in terms of 
cellular phenotypes that are derived primarily from image analysis, but is not so 
limited. Phenotypic stimulus response curves of this invention may contain data 
obtained primarily from non-image sources. Of course, one of ordinary skill in the art 
would recognize other variations, modifications, and alternatives. 
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CLAIMS 

what is claimed is: 

5 1, A method for determining the response of cells to multiple levels of a 
stimulus, the method comprising: 

for each of the multiple levels of stimulus, obtaining a pluraUty of feature 
values, at least some of which characterize the phenotype of cells exposed to the 
particular level of the stimulus, to thereby produce a separate quantitative phenotype 
10 of the cells at each level of stimulus; and 

identifying a path through the separate quantitative phenotypes of cells 
exposed to the stimulus, 

2. The method of claim of 1, wherein the stimulus is selected from the group 
15 consisting of exposure to a chemical compound, exposure to a biological agent, 
exposure to electromagnetic radiation, exposure to particle radiation, exposure to an 
electrical or magnetic field or force, exposure to a mechanical field or force, and 
combinations thereof. 

20 3, The method of claim of 1, wherein the stimulus is exposure to a chemical 
agent. 

4. The method of claim of 1 , fiirther comprising comparing the path to a different 
path produced for a different stimulus to which cells were exposed at multiple levels. 

25 

5. The method of claim of 1, wherein at least some of the feature values are 
obtained from an image of the cells. 

6. The method of claim of 1, wherein at least one of the feature values 
30 characterizes cell morphology. 

7. The method of claim of 1, wherein at least one of the feature values 
characterizes a statistical property of the cells. 

35 8. The method of claim of 1, wherein at least one of the feature values is a 
biological classification of the cells. 
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9. The method of claim of 8, wherein the biological classification specifies a cell 
cycle state of the cells. 

10. The method of claim of 1, wherein the multiple levels of sthnulus are multiple 
5 concentrations of a therapeutic of potential therapeutic. 

11. The method of claim of 1, further comprising presenting a graphical 
representation of the path. 

10 12. The method of claim of 11, wherein the graphical representation is provided 
along one or more principle components obtained via a principle component analysis. 

13. The method of claim of 1, further comprising providing quantitative 
phenotypes for each of multiple stimuli; and 
15 using all the quantitative phenotypes to provide a reduced-dimensionality 

space in which to depict the path. 

14/ The method of claim 1, wherein the multiple levels of stimulus are multiple 
times after an initial exposure to the stimulus. 

20 

15. A computer program product comprising a machine readable medium on 
which is provided program instructions for determining the response of cells to 
multiple levels of a stimulus, the program instructions comprising: 

program code for obtaining, for each of the multiple levels of stimulus, a 
25 plurality of feature values, at least some of which characterize the phenotype of cells 
exposed to the particular level of the stimulus, to thereby produce a separate 
quantitative phenotype of the cells at each level of stimulus; and 

program code for identifying a path through the separate quantitative 
phenotypes of cells exposed to the stimulus. 

30 

16. The computer program product of claun of 1 5, wherein the stimulus is selected 
from the group consisting of exposure to a chemical compound, exposure to a 
biological agent, exposure to electromagnetic radiation, exposure to particle radiation, 
exposure to an electrical or magnetic field or force, exposure to a mechanical field or 

35 force, and combinations thereof. 
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17, The computer program product of claim of 15, further comprising program 
instructions for comparing the path to a different path produced for a different 
stimulus to which cells were exposed at multiple levels. 

5 18. The computer program product of claim of 15, wherein at least some of the 
feature values are obtained from an image of the cells. 

19. The computer program product of claim of 15, wherein at least one of the 
feature values characterizes cell morphology. 

10 

20. The computer program product of claim of 15, wherein at least one of the 
feature values characterizes a statistical property of the cells. 

21. The computer program product of claim of 15, wherein at least one of the 
1 5 feature values is a biological classification of the cells. 

^ 22. ' The computer program product of claim of 15, further comprising program 
code for presenting a graphical representation of the path. 

20 23. The computer program product of claim of 15, further comprising: 

program code for providing quantitative phenotypes for each of multiple 
stimuli; and 

program code for using all the quantitative phenotypes to provide a reduced- 
dimensionality space in which to depict the path. 

25 

24. The computer program product of claim 15, wherein the multiple levels of 
stimulus are multiple times after an initial exposure to the stimulus. 

25. An apparatus for determining the response of cells to multiple levels of a 
30 stimulus from images of the cells, the apparatus comprising: 

an interface configured to receive the images of the cells that have been 

■ 

exposed to said multiple levels of a stimulus; 

a memory for storing, at least temporarily, some or all of the images; and 
one or more processors in conmumication with the memory and designed or 
35 configured to (i) obtain from said images a plurality of feature values, at least some of 
which characterize the phenotype of cells exposed to the particular level of the 
stimulus, to thereby produce a separate quantitative phenotype of the cells at each 
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level of stimulus; and (ii) identify a path through the separate quantitative phenotypes 
of cells exposed to the stimulus. 

26. The apparatus of claim of 25, wherein the stimulus is selected from the group 
5 consisting of exposure to a chemical compound, exposure to a biological agent, 
exposure to electromagnetic radiation, exposure to particle radiation, exposure to an 
electrical or magnetic field or force, exposure to a mechanical field or force, and 
combinations thereof. 

10 27. The apparatus of claim of 25, wherein the one or more processors are further 
designed or configured to compare the path to a different path produced for a different 
stimulus to which cells were exposed at multiple levels. 

28. The apparatus of claim of 25, whereiu the feature values characterizes one or 
15 more of cell morphology, a statistical property of the cells, and a biological 

classification of the cells. 

29. The apparatus of claim of 25, further comprising a display for presenting a 
graphical representation of the path provided by said one or more processors. 

20 

30. The apparatus of claim of 29, wherein the graphical representation is provided 
along one or more principle components obtained via a principle component analysis. 

31. A method for determining whether a first compound and a second compound 
25 act on cells by a related mechanism of action, the method comprising: 

for each of multiple concentrations of the first compound, obtaining a plurality 
of feature values characterizing the phenotype of cells exposed to the particular 
concentration of the first compound, to thereby produce a plurality of first 
concentration-specific phenotypes; 
30 identifying a first path through the first concentration-specific phenotypes of 

cells exposed to the first compound; 

for each of multiple concentrations of the second compound, obtaining a 
plurality of feature values characterizing the phenotype of cells exposed to the 
particular concentration of the second compoimd, to thereby produce a pluraUty of 
35 second concentration-specific phenotypes; 

identifying a second path through second concentration-specific phenotypes of 
cells exposed to the second compound; and 
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comparing the first and second paths, wherein a degree of similarity between 
the paths corresponds to a degree of similarity in the mechanism of action of the first 
and second compoimds. 

5 32. The method of claim 31, wherein at least one of the first and second 
compoimds is a known therapeutic or potential therapeutic. 

33. The method of claim 31, wherein the multiple concentrations of the first 
compoxmd vary from lowest to highest by a factor of at least about two. 

34. The method of claim 31, wherein the multiple concentrations of the first 
compound include at least five separate concentrations of the first compound. 

35. The method of claim 31, wherein the multiple concentrations of the first 
15 compound include at least eight separate concentrations of the first compound. 

36. The method of claim 31, wherein obtaining a plurality of feature values 
characterizing the phenotype of cells exposed to the particular concentration of the 
first compound comprises analyzing images of a population of cells exposed to the 

20 particular concentration of the first compound. 

37. The method of claim 31, wherein the pluraUty of feature values include 
numeric values characterizing one or more of the following cellular components: 
DNA, Golgi, cytoskeletal components, and the plasma membrane. 

25 

38. The method of claim 31, wherein the plurality of feature values include 
numeric values characterizing one or more of the following cellxilar components: 
DNA, Golgi, and tubulin. 

30 39. The method of claim 3 1 , wherein identifying the first path comprises analyzing 
the first concentration-specific phenotypes via oiie or more of the following 
techniques: principal component analysis, linear and non-linear discriminant analysis, 
multidimensional scaling, and projection pursuit techniques. 

35 40. The method of claim 3 1 , wherein identifying the first path comprises analyzing 
the first concentration-specific phenotypes using principal component analysis. 
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41. The method of claim 31, wherein comparing the first and second paths 
comprises graphically depicting the first and second paths together. 

42. The method of claim 41, wherein the graphical depiction presents the first and 
5 second paths in a space defined by principal components. 

43. A computer program product comprising a machine readable medium on 
which is provided program instructions for deteiinining whether a first compound and 
a second compound act on cells by a related mechanism of action, the program 

10 instmctions comprising: 

program code for obtaining, for each of multiple concentrations of the first 
compound, a plurality of feature values characterizing the phenotype of ceUs exposed 
to the particular concentmtion of the first compound, to thereby produce a plurality of 
first concentration-specific phenotypes; 
15 program code for identifying a first path through the first concentration- 

specific phenotypes of cells exposed to the first compound; 

program code for obtaining, for each of multiple concentrations of the second 
compound, a plurality of feature values characterizing the phenotype of cells exposed 
to the particular concentration of the second compound, to thereby produce a plurality 
20 of second concentration-specific phenotypes; 

program code for identifying a second path through second concentration- 
specific phenotypes of cells exposed to the second compound; and 

program code for comparing the first and second paths, wherein a degree of 
similarity between the paths corresponds to a degree of similarity in the mechanism of . 
25 action of the first and second compounds. 

44. The computer program product of claim 43, wherein at least one of the first 
and second compoxmds is a known therapeutic or potential therapeutic. 

30 45. The computer program product of claim 43, wherein the program code for 
obtaining a plurality of feature values characterizing the phenotype of cells exposed to 
the particular concentration of the first compound comprises program code for 
analyzing images of a population of cells exposed to the particular concentration of 
the first compoimd. 

35 

46. The computer program product of claim 43, wherein the plurality of feature 
values include numeric values characterizing one or more of the following cellular 
components: DNA, Golgi, cytoskeletal components, and the plasma membrane. 
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47. The computer program product of claim 43, wherein the program code for 
identifying the first path comprises program code for analyzing the first concentration- 
specific phenotypes via one or more of the following techniques: principal component 

5 analysis, linear and non-linear discriminant analysis, multidimensional scaling, and 
projection pursuit techniques. 

48. The computer program product of claim 43, wherein the program code for 
comparing the first and second paths comprises program code for gr^hically 

1 0 depicting the first and second paths together. 

49. The computer program product of claim 48, wherein the gr^hical depiction 
presents the first and second paths in a space defined by principal components. 
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Legend 

Famesyl transferase (green), 
Geranylgeranyltransf erase (blue) and 
Mitochondrial inhibitors (pink and gray) 
(See legend for Figure 3A for additional details) 
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Legend 

Actin inhibitors (red) 

Tubulin dq>ol)maerizers (purple) 

(See legend for Figure 3A for additional details) 
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Legend 

ER-Ca+2 ATPase inhibitors (black) 

G protein activators (blue) 

Kinase inhibitors (brown and gray) 

(See legend for Figure 3 A for additional details) 
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Legend 

Red data points are from one run 
Blue from a run 3 weeks later 
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Compoimds that have biochemical activity against Cytokinetics' target 
Legend 

Primary hits (red) 
Optimized hits (blue) 
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Legend 

Zoomed-in figure of the compoimds in PC A space in Figure 1 1 . 
Primary hits (red) 
Secondary analogs (black) 
Optimized hits (blue) 
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